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Abstract 

General properties of quantum systems which interact with stochastic environ- 
ment are studied with a strong emphasis on the role of physical symmetries. The 
similarity between the fidelity which is used to characterize the stability of such a 
systems and the Wilson loop in QCD is demonstrated, and the fidelity decay rates 
are derived. The consequences of existence of the symmetry group on the statis- 
tical properties of the system are analyzed for various physical systems - a simple 
quantum mechanical system, holonomic quantum computer and Yang-Mills fields. 



Introduction 

Transition between classical and quantum behaviour is an important phenomenon in 
different fields of physics. Sometimes this transition is associated with the processes of 
decoherence which are commonly described as interaction between quantum system and 
environment in the state of thermodynamical equilibrium [HI E3 \M\ ■ It is known that the 
environment can be described as the stochastic ensemble of classical fields El U21 El El- 
Therefore it is important to analyze the behavior of quantum systems which interact with 
random fields and to relate the statistical properties of random fields to the physical 
properties of the environment, for example, the properties of symmetry. Implementations 
of the model of stochastic environment with definite symmetry properties are very diverse - 
for example, the stochastic vacuum of Yang-Mills field which is known to possess confining 
properties, or the errors in quantum computations. The stability of quantum motion is 
usually characterized in terms of the fidelity. In this paper the expression for the fidelity 
of holonomic quantum computations is derived, and the analogy between the theory of 
holonomic quantum computations and the theory of gauge fields is demonstrated. It is 
shown that the roots of this analogy lie in the existence of invariance group of the system. 
This analogy is used to introduce the concept of the fidelity in QCD. It is shown that the 
fidelity decay rate essentially depends on the " string tension" . Invariance of the stochastic 
ensemble of classical fields under the action of the symmetry group is used to simplify the 
calculations. 



1. Multiplicative integrals with random variables 

The quantities which are analyzed in this paper are most naturally expressed in terms 
of the expectation values of multiplicative integrals over some subspaces in the config- 
uration space of the problem. Multiplicative integral of operator field A (x) over some 
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manifold D of dimensionality d is defined as pQ: 
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where P is the path-ordering operator, I is the identity operator, dV is the volume ele- 
ment on D and AVk is the small finite element of volume on D. Both dV and A (x) can 
be tensors, but their product should be a scalar quantity. 

Generally there is no expression in closed form for the expectation value of multi- 
plicative integral of an arbitrary random field A (x) , but a good approximation can be 
found when the size of the manifold D considerably exceeds the correlation length / corr of 
random field. In this case the manifold D can be split into a large number m of subman- 
ifolds Dk, k = 1 . . . m with the sizes still much greater than the correlation length. The 
multiplicative integral (0) is then rewritten as the product of multiplicative integrals over 
each submanifold Dk. As the sizes of submanifolds Dk exceed the correlation length of the 
random field, multiplicative integrals over them can be treated as statistically indepen- 
dent variables, and expectation value of their product can be calculated as the product 
of expectation values of each factor. Each multiplicative integral over submanifold Dk is 
then expanded in powers of A (x). The result of these transformations is: 
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where we have used the fact that the operation of path-ordering is linear and can be 

therefore exchanged with averaging. The function C (x, y) = A (x) A (y) — A (x) A (y) is 
the second-order correlation function of random field A (x) . Here we will consider only 
statistically homogeneous random fields. A random field is statistically homogeneous if 
its correlation functions depend only on the differences of their arguments 0. In this 
paper the following form of the correlation function is assumed: 

(7(x,y)=C7(x-y), / (0) = 1 (3) 

The function / (x) is significantly large only when |x| < / corr . The second-order term in 
the expansion (J2J) is then approximated to the first order in the volume AVk — J dV of 

D k 

submanifold D k as P f f dVidA^A (xi) A (x 2 ) ~ l/2l^ OII CAVk- High-order terms in this 

D k D k 

approximation are estimated as P J ... J dV 1 ... dV { A (x x ) ...A (x,) ~ C AV k . 

D k D k 

High-order terms can be neglected if the following condition holds: 

C' 1/2 C r «l (4) 

This condition is called the condition of gaussian dominance El- For the sake of 
simplicity we will assume that this condition holds. The contribution of high-order terms 
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can be taken into account without any difficulties. Under these assumptions the expression 
(J2J) can be rewritten in a compact form using the definition of multiplicative integral 
and the fact that the number m is large: 



U=\{(i + iA (x fc ) AV k - 1/2 l d COTT CAV k 
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(5) 

If the stochastic ensemble of fields A (x) is invariant under global transformations T 
which build an irreducible representation of some group G, the expectation value of the 
multiplicative integral (JTJ is proportional to the identity. Invariance of the stochastic 
ensemble is understood as the following identity: 
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where DP 
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is the probability measure on the stochastic ensemble of fields A (x). 
This assertion is proved by applying an arbitrary transformation from the group G to the 
expectation value of the product of an arbitrary number of field variables A (x): 
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Thus the expectation value of the product of an arbitrary number of field variables com- 
mutes with an arbitrary transformation from the irreducible representation of the group 
G. But then the expectation value of the multiplicative integral also commutes with 
all such transformations. This is simply proved by expanding the multiplicative integral 
in powers of its argument and applying the transformation to each summand. Schur's 
lemma in group theory states that if the linear operator commutes with all transforma- 
tions from the irreducible representation of the group it is proportional to the identity [3]. 
The assertion is proved. It is evident that in this case path-ordering is no longer necessary 
in ©. 



2. The fidelity. A simple example: quantum system with random hamil- 
tonian 

The fidelity of an arbitrary quantum system is defined as the scalar product of the 
state vectors which correspond to perturbed and unperturbed states of the system. If 
there are random variables in the problem, the expectation value of the square of the 
absolute value of the fidelity is considered [TTJ H2J [TH] : 

A = <0o|0i>, h = W = Tr (p 1 /) ) (8) 

where \<f>i), \4>o) are the state vectors of perturbed and unperturbed systems and p 1 and 
p° are the density matrices of perturbed and unperturbed systems. The fidelity is used to 
characterize the stability of quantum systems with respect to perturbations [TT | IT2" | IT3]. 
In this section a simple example of quantum system which interacts with stochastic 
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process with definite symmetry properties is considered. On this simple example we will 
demonstrate that it is possible to use the expectation value of the fidelity fi to estimate 
the fidelity decay rate and that decay rates of the expectation value of the fidelity j\ and 
the fidelity f'2 are equal up to some coefficient which is close to one. 

Statistical averaging in quantum mechanics is performed on the density matrix of the 
system. Evolution of the density matrix is described by the Liouville equation (h = 1): 

dtPij = i (pikHkj — Hikpkj) (9) 

We will suppose that the hamiltonian Hij is the sum of the unperturbed stationary hamil- 
tonian and the random perturbation, which is the stationary stochastic process [T3*t ITlj . 
The total number of independent state vectors of the system is assumed to be N. We 
will also suppose that the matrices of the perturbation part of the hamiltonian form the 
gaussian orthogonal ensemble [2 which is invariant under the action of the orthogonal 
group O(N): 

H tj = H% + SHij, H% = diag (E { ) , 

6Hij 5H kl (t 2 ) = o 2 (b ik b jX + 5 u 5 jk ) fih-h), Jdtf (t) = r corr (10) 



where r corr is the correlation time of the perturbation. Solution of the Liouville equation 
(J5jl can be represented in terms of multiplicative integral after introducing the linear 
operator Xijki = Hj5ik — H ik 5ij on the space of [N x N] matrices. The Liouville equation 
is rewritten as d t pij = iXijkiPki, an d its formal solution is obtained by applying the 
multiplicative integral of the operator X^i to the density matrix of the initial state: 



p(t) = Pexp(i jdZX(Z)\p(p) (11) 



Now it is possible to use the results of section 1 and to calculate the density matrix which is 
averaged over all implementations of stochastic process SHij. The manifold D is the time 
interval from to t, and the condition of gaussian dominance (0J) is crr corr 1. The expec- 
tation value of the operator X^i is X^ k i = H^8 ik — H® k 5ij, the second-order correlation 
function is Cy OTn = X^X klmn — Xij k iX klmn = 2cr 2 (6j n 5 im (N + 1) — 5 ni 8j m — 5 nm 5ij). 
The density matrix is obtained as the expectation value of the density matrix in (jllj) 
according to (J3J): 

p(t)=P exp [iXt - 1/2 r^Ctj p (0) (12) 

Note that the operator X is invariant under the transformations from the orthogonal group 
O (N) , but these transformations do not build irreducible representation, and the second- 
order correlation function is not proportional to the identity. However the expression (|T2*|) 
can be calculated in closed form and the final result is: 

PH (*) = (PH (°) - 1 / N <y ex P (« ( e j -Ei)t- a 2 Nr covv t) + 1/N 5 i3 (13) 

The density matrix of the unperturbed system is if) = pij (0) exp [i (Ei — Ej) t). The 
fidelity is f% (t) = p^ if) p^ [t) = (1 — 1/N) exp (— <7 2 iVT corr t)-|-l/iV. The expectation value 
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of the fidelity f\ is obtained from (jSJ) in a similar way using the second-order correlation 
function of the hamiltonian H^H^j — Hn-H^j = a 2 (N + 1) <5y, and the result is fi (t) = 
exp (—1/2 <7 2 T corr (N + l)t), which is in a good agreement with the exact expression for 
f2 (t). Thus the fidelity decay rate which is obtained in exact but somewhat complicated 
way and those which is obtained by simple averaging differ only by the factor N/ (N + 1), 
which can be explained by randomly changing phase of the fidelity f\. The asymptotic 
behaviour is different, but it is not considered here. Further in this paper the expectation 
value of the fidelity fi is used to estimate the fidelity decay rate. 



3. Wilson loop 

Wilson loop is commonly used to test the confining properties of Yang-Mills fields. 
Wilson loop is defined as the trace of the multiplicative integral of gauge field A^ (x) over 
some closed contour 7 in Minkowski space E] : 



W(i) = l/iV c TrPexp ig / dx^A 



(14) 



where the hat symbol denotes operators in the colour space, N c is the number of colours 
and g is the coupling constant. If the "area law" holds for the Wilson loop W (7), that is 
W (7) ~ exp (—aS 7 ) , 5 7 being the minimal area of the surface spanned over the contour 
7, quarks are said to be tied by a string with the constant "tension" a @||. If the 
gauge fields A^ (x) form stochastic ensemble, the expectation value of the Wilson loop 
is considered jSJEl- In the case of stochastic vacuum in QCD one usually chooses the 

curvature tensor of gauge field F^ v = d ii A v — d v A^ — ig A M , A u as a random variable 

jHUE]. For topologically trivial fields the multiplicative integral over the contour 7 in (fT4*|) 
can be represented as the multiplicative integral over the surface spanned over the contour 
7 by applying the non-Abelian Stokes theorem El E] : 



W{i) = l/N c TrPexp 



ig J dS^F^ 



(15) 



where F^ v (y) = U (x, y) F M „ (y) U (y, x) is the shifted curvature tensor, 
U (x, y) = Pexp ig J dx^A^ ) and x is an arbitrary point on the contour 7. The di- 
mensionality d of the manifold D in (0) is 2 in this case, and the condition of gaussian 
dominance (j3J) for the random field F^ v (x) is: 
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V < 1 



(16) 



QCD is by definition invariant under local gauge transformations. The transformation 
law for the curvature tensor F^ is F^ v (x) — > T\s)F^ (x) T(x), where T (x) e SU (3). As 
the gauge symmetry should be also preserved for the random field F^ v (x), the probability 
measure on F^ v (x) should be invariant under gauge transformations. Now the results of 
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section 1 can be applied to calculate the expectation value of the Wilson loop. The 
random field is the field of curvature tensor, the manifold D is the surface spanned over 
the contour 7, dV is the element of area, the invariance group G is the group of local gauge 
transformations, which obviously includes the subgroup of global gauge transformations. 
The final expression for the expectation value of the Wilson loop according to (J5J) is: 

WW) = exp (-1/2 g 2 Tnl mx S y ) (17) 

Thus the Wilson area law for the Wilson loop holds, and the gaussian-dominated stochas- 
tic vacuum possesses confining properties. The string tension a is also obtained from (JT7J): 
a = 1/2 g 2 l 2 COII F 2 . 

4. Holonomic quantum computations 

Holonomic quantum computations are implemented as unitary transformations of 
the vector space spanned on eigenvectors corresponding to the degenerate energy level 
of the hamiltonian H which depends on the set of external parameters A M : H (A) = 
U (A) HqW (A), where Z7 (A) is unitary. Typically degeneracy of the energy levels is 
the consequence of existence of the invariance group G of the hamiltonian H (A M ) 10 . 
Wigner's theorem states that if the hamiltonian of the system is invariant under transfor- 
mations T which build a representation of the group G on the Hilbert space of the system, 
that is T~ X HT = H, each energy level is degenerate, and transformations of the vector 
subspaces of the Hilbert space of the system which correspond to the degenerate energy 
levels build irreducible representations of the group G [3]. Only one degenerate energy 
level with eigenvectors which build a basis of an irreducible representation is considered 
further. Holonomic quantum computations are performed by changing the parameters 
A^ adiabatically [TO]- In this case dynamical effects can be neglected, and an arbitrary 
quantum gate is then described by the multiplicative integral over some contour 7 in the 
space of external parameters: 



f 7 = Pexp [i / dAM„(A) 




where A W y (A) = (0j (A)| A M (A) \<pj (A)) = — i ((pi (A)| \<pj (A)) is the adiabatic connec- 
tion on the space of external parameters and |0j(A)), i = 1...N are the eigenvectors 
which correspond to the degenerate energy level. 

It is interesting to note the close analogy between the theory of gauge fields and the 
theory of holonomic quantum computations. For example, in the theory of gauge fields 
gauge transformations correspond to local transformations of basis vectors in the colour 
space, and the transformation law for the gauge field is (x) — > T _1 (x) A M (x) T (x) — 
iT^ 1 (x) d^T (x). Physical observables obviously do not depend on the choice of basis vec- 
tors, and the consequence is the gauge invariance of the theory. For holonomic quantum 
computations it is possible to consider unitary transformations of basis vectors which cor- 
respond to the degenerate energy level. These transformations in general depend on the 
external parameters A M . The adiabatic connection (A) transforms under local trans- 
formations of basis vectors as A^y (A) — > —i (0j (A)| T _1 (A) gfjrT (A) \<pj (A)) = 
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N 

-i £ (^(A)|T- 1 (A)|0 fc (A))(^(A)|^r|^(A))(^(A)|T- 1 (A)|^(A)) = 
fc,/=i 

Tg 1 (A) ,H (A) Tij (A) - (A) ^T kj (A), or, equivalently, 

(A) -> f - 1 (A) i„ (A) f (A) -«f - 1 (A) (A) which is evidently the same as the trans- 
formation law for the gauge field. The principle of gauge invariance can be reformulated in 
this case as invariance of physical observables with respect to the choice of irreducible rep- 
resentation of the group G among all equivalent irreducible representations. Transitions 
between equivalent irreducible representations correspond to unitary transformations of 
the basis vectors of the representations [4^ . 

The fidelity of holonomic quantum computations with respect to errors in control pa- 
rameters is defined as the scalar product of the state vectors which correspond to desired 
and actual operation. Using the general expression for the quantum gate (|T£|) . the fidelity 
of the holonomic quantum computation can be expressed in terms of the multiplicative 
integral [TT|IT2]: 

/=(V>o|Pexp \i I dAMj |Vo) (19) 




57 

where £7 is the contour obtained by traveling forward in the path 71 which corresponds to 
desired operation and backward in the path 7 2 which corresponds to actual operation and 
|t/>o) is the initial state of the quantum register. As the fidelity is usually close to unity 
(otherwise holonomic quantum computations will lead to unpredictable results and will 
be therefore useless), the distance between the contours 71 and 72 should be very small. In 
this case it is possible to rewrite (fT9*j) in the form where integration is performed over only 
one contour. This is achieved by continuously deforming the contour 7 2 to the contour 
71 and simultaneously changing the vector on it in such a way that the multiplicative 
integral in ([T9*|) is not changed: 

Pexp j i J dVA^ J = Pexpli J (i M + 5A^j J (20) 
\ 72 / \ 71 / 

To obtain the correcting term SA^ (X u ) we multiply the equation (j20j) from the left by 
P exp I J dX^A^ ) ) . The left side is transformed by applying the non-Abelian Stokes 

\7l 

theorem: 



(21) 




Pexp i J dX^A^ Pexp i J dX^A^ 

\ 71 / / \ 72 

Pexp j i j dX^A^ J = Pexp ( i j dX»F^5X v (A) 

V $1 J V 71 

where S X v is the distance between corresponding points on the contours 71 and 72 and 
is the shifted curvature tensor defined as in section 3 with the gauge field A M (x) replaced 
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by the adiabatic connection A M (A). The right side of (j2*U|) is transformed as: 



-1 



Pexp J dX^A^j Pexp J dX^ 1 (A^ + SA^j 



(22) 



Pexp [i J dX^U- 1 (A) 5AJJ (A) 



Ti 



where U (A) = Pexp ^ J dX^A^j is the multiplicative integral of the adiabatic connection 

A^ (A) over the piece of the contour which is limited by the initial point and the point 
A. The correcting term SA^ in ()2U|) is then obtained by comparing the arguments of the 
multiplicative integrals in (}22|) and (|2~T|) : 



5A tl {X) = F fMV (X)5X 1 ' (23) 

The expression (I19j) for the fidelity of holonomic quantum computations is then rewritten 
as: 

/ = <-0 o | Pexp ^ ^ dX^U- 1 (A) (A) 17 (A) j |^ > (24) 

The error SX is usually random. It is therefore important to predict the fidelity decay 
rate from the statistical properties of the errors. As in section 2, the expectation value 
of the fidelity will be used to estimate the fidelity decay rate. The expectation value 
of the multiplicative integral in can be calculated using the expression Q. The 
manifold D is now the contour 71, and dV is the element of length on it. One can expect 
that the probability measure on 5A^ (A) does not depent on the choice of the irreducible 
representation, hence the invariance group is G, and the multiplicative integral in ()19|) 
according to section 1 is proportional to the identity. Therefore the fidelity does not 
depend on the initial state \ipo), and the final expression for the expectation value of the 
fidelity is: 



/ = exp -1/2 A corr / 5A^5A v n^dX l 



71 



(25) 



exp ( -1/2 A corr / dX u F xa (A) F vj3 (A) ^A^A^n* 
71 



where n x is the directing vector of the path 71 and A cor r is the correlation length of the 
error ^A^ (A) which indicates the relative frequency of errors along the contour 71. The 
expression (|23|) shows that the stability of quantum computation essentially depends on: 



the duration and the closeness of errors (the factor A corr in (J2SJ)) 
the curvature tensor on the space of external parameters F^ u 
the expectation value of the square of the magnitude of errors 
the length of the contour 7! in the space of external parameters 
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5. Fidelity in QCD 

Using the analogy between the theory of gauge fields and the theory of holonomic 
quantum computations which was pointed out in section 4 we define the fidelity of quark 
as the scalar product of state vectors in the colour space \qi) and \q%) which correspond 
to perturbed and unperturbed motion: 

/=<?i|?2>, ki> , |g 2 ) G C 3 , ( qi \ qi ) = (q 2 \q 2 ) = 1 (26) 

It is now possible to apply the results which were obtained in previous sections and to 
obtain the fidelity decay rates for different types of quark motion. First we consider the 
motion of a coloured quark in different paths 71 and 72 in the QCD stochastic vacuum. 
The paths start from the point x and join in the point y. In the initial point x the state 
vector is \qo). I n the limit of very massive quarks E] evolution of the state vectors in 
the colour space is described by the multiplicative integral introduced in ([14)1: 

\q k ) = Pexp\ig [ dx^A^ \q ) = U (7*) \q ) , k = 1, 2 (27) 



7/b 



The operators U (71) and U (72) are unitary because is hermitan. Taking this into 
account, one can rewrite the expression for the fidelity as the multiplicative integral over 
the closed contour 7 = 7172 obtained from the paths 71 and 72 by travelling from the 
point x to the point y in the path 71 and back to the point x in the path 72 : 



f=(qo\U (71) ■ U+ (72) \qo) = (qo\U (7) |g ) (28) 

As in the section 2 we use the expectation value of the fidelity to estimate the fidelity 
decay rate. All calculations which are necessary have been already made in section 3. The 
multiplicative integral in (|2*8j) is proportional to the identity due to the colour neutrality 
of the stochastic vacuum. The expression for the fidelity then does not depend on the 
initial state \qo) because of the normalization. The condition of gaussian dominance is 
([16]) . The final expression for the fidelity of the quark moving in the gaussian-dominated 

stochastic vacuum is: 

/ = exp(-l/2^ 2 / c 2 orr F 2 5 7 ) (29) 

where S" 7 is the area of the surface spanned over the contour 7. Thus for the gaussian- 
dominated stochastic vacuum the fidelity for the quark moving in different paths decays 
exponentially with the area of the surface spanned over the paths, the decay rate being 
equal to the string "tension" a. This hints at the close connection between the stability 
of quark motion and quark confinement. 

Another possible situation, which is more close to the standard treatment of the fi- 
delity, is realized when 7! and 72 are two random paths in Minkowski space which are very 
close to each other. The corresponding expression for the fidelity is similar to ([28]) . but 
now the averaging is performed with respect to all random paths which are close enough. 
All the necessary calculations are similar to those in the section 4. The invariance group 
is again the group of local gauge transformations. The final expression for the fidelity in 
this case is: 

/ = exp ^-1/2 g 2 l COTT J dx u F xa F vP v*J^5^ (30) 
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where 5x a is the deviation of the path 72 from the path 71, v x is the four-dimensional 
velocity and / corr is the correlation length of perturbation of the quark path expressed in 
units of world line length. For example, if unperturbed path is parallel to the time axis 
in Minkowski space, the quark moves randomly around some point in three-dimensional 
space. The fidelity in this case decays exponentially with time, as would be expected. 

Conclusions 

In this paper the interactions between different quantum systems and the stochastic 
environment was studied. It was shown that the approach described in section 1 is uni- 
versal and can be applied to a wide class of systems in different branches of physics - from 
simple quantum systems to quantum computer and gauge field theories. It is shown that 
existence of the invariance group of the environment to a large extended determines the 
behavior of quantum system. This phenomenon is intensively investigated in the Random 
Matrix Theory [2|, where the statistical properties of eigenvalues of operators are consid- 
ered. The stochastic ensemble of fields which is invariant under the action of some group 
can be regarded as the straightforward generalization of orthogonal or unitary ensembles 
(GOE, GUE) in the Random Matrix Theory. However it was sufficient to require in- 
variance under the global transformations only to obtain such an important properties of 
stochastic ensemble as proportionality of the correlators to the identity. The consequences 
of invariance under local gauge transformations, as well as the role of gauge invariance 
in the emergence of chaotic behaviour require further investigations. It was also shown 
that the fidelity decay rate for quarks moving in QCD stochastic vacuum is equal to the 
" string tension" which characterizes the confining properties of Yang-Mills fields. There- 
fore a relation exists between the property of quark confinement in stochastic vacuum and 
the stability of quark motion. This fact can be related to possible mechanisms of quark 
confinement. 
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